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In this lecture we give a brief review of chiral Random Matrix Theory (chRMT) and its 
applications to QCD at nonzero chemical potential. We present both analytical arguments 
involving chiral perturbation theory and numerical evidence from lattice QCD simulations 
showing that correlations of the smallest Dirac eigenvalues are described by chRMT. We 
discuss the range of validity of chRMT and emphasize the importance of universality. For 
chRMT's at /i 7^ we identify universal properties of the Dirac eigenvalues and study the 
effect of quenching on the distribution of Yang-Lee zeros. 

1. INTRODUCTION 

There are strong indications from hadronic phenomenology that chiral symmetry is 
broken spontaneously at low temperatures. However, mainly through lattice QCD simu- 
lations, it has become generally accepted that chiral symmetry is restored above a critical 
temperature of about 140 MeV (see reviews by DeTar Ukawa 0, Smilga [|3] and 
Karsch || for recent results on this topic). The situation at nonzero chemical potential 
is much less clear. Because of the phase of the fermion determinant it is not possible to 
perform Monte-Carlo simulations, and it seems that one has to give up calculations based 
on first principles. It is certainly not possible to make progress by a gradual improvement 
of existing techniques. Instead, radically different approaches have to be developed. 

The main source of the problem is the loss of Hermiticity at nonzero chemical potential. 
In this lecture we study this problem by means of a much simpler chiral Random Matrix 
Theory (chRMT) with the global symmetries of the Dirac operator [H-H. This model 
has the remarkable property that, in spite of the fact that it can be solved analytically, it 
cannot be solved numerically. In the broken phase, numerical convergence is only obtained 
for an exponentially large ensemble ||. If one is not able to develop an algorithm for this 
simple model, progress in QCD will be elusive. 

Originally, chRMT was introduced to describe the correlations of QCD Dirac eigen- 
values. Our interest in the Dirac spectrum is based on the relation between the chiral 
condensate, S, and the spectral density per unit of space time volume, V, given by 
£ = Iim7rp(0)/U ||10|| . Note that the thermodynamic limit has to be taken before the 
chiral limit. The spectral density is defined by p(\) = J2k${^ ~ ^k), where the are the 
eigenvalues of the Dirac operator. Because of this relation the eigenvalues near zero are 
spaced as l/p(0) = ir/YiV. In order to study the approach to the thermodynamic limit it 
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is natural to introduce the microscopic limit in which u 
and the microscopic spectral density 
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XVTj is kept fixed for V — > oo, 
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Our claim is that the microscopic spectral density of the QCD Dirac operator is given by 
chRMT. This has been confirmed both by analytical arguments using partially quenched 
chiral perturbation theory |IT| , and lattice QCD |I2|-|I6| and instanton liquid |T7| , |T8l 



and lattice QCD JT2|- 

simulations. In addition to this, correlations in the bulk of the spectrum are given by 
chRMT JTSH^JTSI as well. Random matrix theory can describe only those observables 
which are universal, i.e., which are stable against large deformations of the random matrix 
ensemble. By now it has been well established that both the microscopic spectral density 
and the spectral correlations are strongly universal f22|^2~9|| , | 30fl . Before discussing the 
random matrix model at nonzero chemical potential, we wish to establish the domain of 
applicability of chRMT at fi = 0. Starting from partially quenched Chiral Perturbation 



Theory |3T| we will argue that, below a scale of Aq C d/VV, eigenvalue correlations are 
given by chRMT. This scale is the equivalent of the Thouless energy in mesoscopic physics 
l2|-011 for recent reviews). The existence of such a scale has been confirmed by lattice 
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QCD 0|T^] and instanton liquid simulations. The interpretation is that the classical 
motion of the quarks in the background gauge fields is chaotic. 

Before proceeding to the main body of this talk, let me stress that there are two different 
types of applications of RMT. First, as the simplest model of a universality class, and 
second, as a schematic model for a complex system. Typical examples in the first class are 
related to eigenvalue correlations ||32|| . Perhaps, the most famous example in the second 
class is the Anderson model for localization ||35|| . In this lecture we will discuss a schematic 
RMT model for the chiral phase transition at nonzero chemical potential. 

2. CHIRAL RANDOM MATRIX THEORY 



In this section we will introduce an instanton liquid p6| , |37| inspired chiral RMT for the 
QCD partition function. With only its global symmetries as input but otherwise Gaussian 
distributed matrix elements the chRMT for Nf flavors and topological charge v is defined 

by 1® 



m N . 



N f 

f DW J] det(X> + m/)e~^ 



where 
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and W is a n x m matrix with v — \n — m\ and N = n + m. We interpret N as the 
(dimensionless) volume of space time. The matrix elements of W are either real (j3 = 1, 
chiral Gaussian Orthogonal Ensemble (chGOE)), complex (/3 = 2, chiral Gaussian Unitary 
Ensemble (chGUE)), or quaternion real {(3 = 4, chiral Gaussian Symplectic Ensemble 
(chGSE)). As is the case in QCD, we assume that v does not exceed vN, so that, to a 
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good approximation, n = N/2. The Dyson index (3 of a physical system is determined by 
its anti-unitary symmetries. If the square of the anti-unitary symmetry operator is equal 
to the identity, we have (3 = 1. If its square is equal to minus the identity, we have (3 = 4. 
If there no anti- unitary symmetries, the value of (3 = 2. 

In this model chiral symmetry is broken spontaneously with chiral condensate given 
by S = lini7v_ + oo irp(0)/N. For complex matrix elements (J3 = 2), which is appropriate 
for QCD with three or more colors and fundamental fermions, the symmetry breaking 
pattern is SU(N f ) x SU(N f )/SU(N f ) i,fH • 

The average spectral density that follows from (|2]) has the familiar semi-circular shape. 
The microscopic spectral density can be derived from the limit ([!]) of the exact spectral 
density for finite N. For N c = 3, Nf flavors and topological charge v it is given by p0|j38[] 

Ps(z) = ~ (J 2 a (z) - J a+1 (z)J a ^(z)) , (4) 

where a = Nf + \v\. The spectral correlations in the bulk of the spectrum are given by 
the invariant random matrix ensemble with the same value of (3 pT| , [42f . 

The description of an observable in terms of chRMT is based on universality, i.e. on 
the stability against deformations of the random matrix ensemble. At this moment I 



only mention the work of Akemann, Damgaard, Magnea and Nishigaki [p2] who showed 
that the same microscopic spectral density (Q) is obtained for an invariant probability 
distribution defined by an arbitrary polynomial potential. The essence of their proof is a 
remarkable generalization of the identity for the Laguerre polynomials, lim^oo L n (x/n) = 
J (2y / x), to orthogonal polynomials determined by the probability distribution. 

3. DOMAIN OF VALIDITY 

The domain of validity of chRMT is best discussed within the context of partially 
quenched chiral perturbation theory [51, IT]. This is an effective field theory for the low- 



energy limit of a QCD like theory which, in addition to the usual Nf sea quarks, contains 
k valence quarks with mass m v and their super-symmetric partners with the same mass. 
In this framework it is possible to calculate the valence quark mass dependence of the 
chiral condensate which is defined as 

K(m v ) = 1 / dX(p(X))-^^. (5) 
N J X s + mfj 

The spectral density follows from the discontinuity of E„(m„), 

2ir 2-n 

-(p(A)> = Disc|„ A S(m„) = limS(zA + e) - S(zA - e) = — J>(A + A fc )>, (6) 

k 

where the average (■ ■ •) is with respect to the distribution of the eigenvalues. Similarly, 
the two-point spectral correlation function, given by 

(/?(A)p(A/)) = h DisC l-^,"V=a< 5 (iAfc + m^ + nv) ' (?) 

is related to the scalar susceptibility. Both the valence quark mass dependence of the chiral 
condensate and the scalar susceptibility can be calculated from the pqChPT partition 
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function. In both cases we can identify an important scale where the mass of the Goldstone 
modes containing a valence quark is equal to the inverse length of the box ||3,||J]. Using 



the relation M = (m + m')E/F 2 , where F is the pion decay constant, we find from 
ML = 1 that, in terms of the valence quark mass m v , this scale is given by fT8| , |4"5|| 



E, = m . (8) 

For m v <C E c we have shown that the valence quark mass dependence and the scalar 
susceptibility can be obtained from the zero-momentum component of the pqChPT par- 
tition function and are given by the result for chRMT. The conclusion is that, if pqChPT 
describes correctly the low energy limit of QCD, we have shown that the correlations of 
the Dirac eigenvalues close to zero are given by chRMT. 

Such a picture is well-known from mesoscopic physics. In this context E c is defined 
as the inverse tunneling time of an electron through the sample which is given by E c = 
hD/L 2 , where D is the diffusion constant. Another scale that enters in these systems is 
the elastic scattering time r e . Based on these scales one can distinguish three different 
regimes for the energy difference, 5E, that enters in the two-point correlation function 
IE]: i) the ergodic regime for SE <C E c , ii) the diffusive domain for E c -C SE <C fi/r e 



and iii) the ballistic regime for 5E 3> h/r e . On time scales corresponding to the ergodic 
regime an initially localized wave packet covers all of phase space. In this domain the 
eigenvalue correlations are given by RMT. On time scales corresponding to the diffusive 
domain an initially localized wave packet explores only part of the phase space resulting in 
eigenstates with wavefunctions that are localized in different regions. The corresponding 
eigenvalues are show weaker correlations which are no longer given by RMT. For earlier 
applications of localization theory to the chiral phase transition we refer to H7 . 



Based on these ideas we can interpret the Dirac spectrum as the energy levels of a 
system in four Euclidean dimensions and one artificial time dimension. According to 



the Bohigas conjecture [|48| the eigenvalue correlations are given by RMT if and only if 
the corresponding classical motion is chaotic. We thus conclude that the classical time 
evolution of quarks in the Yang-Mills gauge fields is chaotic. 



These ideas have been tested by means of lattice QCD fll4|,|15 | an d instanton liquid fll8 



simulations. It has been found that the eigenvalue correlations are given by chRMT up 
the predicted scale of E c = F 2 /T,L 2 . 

4. LATTICE QCD RESULTS 

In this section we consider correlations of lattice QCD Dirac eigenvalues. For a Dirac 
operator with a Ua(1) symmetry the eigenvalues occur in pairs ±Afe. Therefore we have 
to distinguish two different regions: the region near zero virtuality and the bulk of the 
spectrum. The £7a(1) symmetry is absent for the Hermitean Wilson Dirac operator. 

The class of random matrix ensembles is determined by anti-unitary symmetries. In the 
case of an SU (2) color group, the anti-unitary symmetries of the Kogut-Susskind (KS) 



and the Wilson Dirac operator are given by ||4"9",l9|, 

[D ks ,t 2 K] = 0, and [ l5 D w , j 5 CKr 2 ] = 0. (9) 
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Because (T2K) 2 = — 1 and (75CK72) 2 = 1, we find a Dyson index (3 = 1 for Wilson 
fermions and (3 = 4 for KS fermions. The KS Dirac operator for two colors is thus 
described by the chGSE. In the absence of a Ua(1) symmetry, the Wilson Dirac operator 
for two colors is described by the GOE ||19||. For three or more colors there are no anti- 



unitary symmetries and the relevant ensembles are the chGUE and GUE for KS and 
Wilson fermions, respectively. 

In order to separate the correlations from the average spectral density, the eigenvalues 
are rescaled according to the local average level spacing. Under the assumption of spectral 
ergodicity, eigenvalue correlations in the bulk of the spectrum can be calculated by spectral 
averaging instead of ensemble averaging. For subtle differences between the two we refer 
to ||15|| . Eigenvalues calculated by Kalkreuter ]50| for the N c = 2 KS and Wilson Dirac 
operators for lattices as large as 12 4 have been analyzed by means of spectral averaging 
ITOipDI. For different values of (3, ranging from strong coupling to weak coupling, it was 
found that correlations are in complete agreement with chRMT for distances as large as 
100 level spacings. Recently, correlations in the bulk of the spectrum were also investigated 



for N c = 3 |2~T| , |5Tf . Agreement with chRMT was even found for values of (3 well above the 
deconfinement transition. For a discussion of global Wilson spectra we refer to p2"||)3"|| . 

Spectral ergodicity cannot be exploited in the study of the microscopic spectral density. 
In order to gather sufficient statistics one has to generate a large number of independent 
QCD Dirac spectra. The analysis of these spectra convincingly shows that the smallest 
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Figure 1. The distribution of the smallest eigenvalue (left) and the microscopic spectral 
density (right) of the Kogut-Susskind Dirac operator for two colors and (3 = 2.0. 



eigenvalues are correlated according to chRMT [p^2| - [l^J16| , p!5[| . Results for 1416 configu- 
rations on 10 4 lattice for (3 = 2 are shown in Fig. (|1]). For more detailed results, 
including results for the two-point correlation function, we refer to the original work. In 
Fig. 4 we show the distribution of the smallest eigenvalue (left) and the microscopic 
spectral density (right). The lattice results are given by the full line. The dashed curves 
represent the random matrix results. We emphasize that the theoretical curves have been 
obtained without any fitting of parameters. The input parameter, the chiral condensate, 
is derived from the same lattice calculations. Recently, the same analysis |L2| was per- 
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formed for (3 = 2.2 and for (3 = 2.5 on a 16 4 lattice. In both greement with the 

random matrix predictions was found |L2] . 

It is an interesting question of how spectral correlations of KS fermions evolve in the 
approach to the continuum limit. Certainly, the Kramers degeneracy of the eigenvalues 
remains. However, since Kogut-Susskind fermions represent 4 degenerate flavors in the 
continuum limit, the Dirac eigenvalues should obtain an additional two-fold degeneracy. 
We are looking forward to more work in this direction. 



5. APPLICATIONS OF chRMT AT p ^ 

In the continuum formulation of QCD the chemical potential enters in the QCD parti- 
tion function by the addition of the term /ryo to the anti-Hermitean Dirac operator, i.e. 
T> — > P + /i7 . In a suitable chiral basis, in which the matrix elements of (0#|7o !</>£) = ^kh 
the corresponding modification in the random matrix partition function (§) is to replace 
the Dirac matrix T> by || 

As is the case in QCD, a nonzero chemical potential violates the ant i- Hermit icity of 
the Dirac operator and its eigenvalues are scattered in the complex plane. For three or 
more colors this results in a complex fermion determinant. The presence of this phase 
makes it virtually impossible to perform Monte-Carlo simulations for the full theory. The 
easy way out would be to ignore the fermion determinant altogether. However, it was 
shown that the quenched approximation fails dramatically p5|-57|. The critical chemical 



potential is determined by the pion mass instead of the nucleon mass. After some earlier 
work by this puzzle was resolved convincingly by Stephanov || within the context 

of chRMT at /i 0. He could show analytically that the quenched approximation is not 
the Nf — > limit of the full theory but rather the Nf —>■ limit of a theory in which the 
fermion determinant is replaced by its absolute value. This theory can be interpreted as 
a theory that contains an equal number of quarks and conjugate quarks. It was shown 
that in this theory chiral symmetry is broken spontaneously by the formation of a quark 
conjugate anti-quark condensate. The critical chemical potential is then determined by 
the mass of the corresponding Goldstone bosons with a net baryon number. 

Because the term fijo does not change the anti-unitary symmetries of the Dirac operator 
the classification for ji = is also valid at fi ^ 0. For example, the Kogut-Susskind Dirac 
operators for N c > 3 and N c = 2 corresponds to the classes (3 = 2 and (3 = 4, respectively, 
whereas the continuum Dirac operator for N c = 2 corresponds to /3 = 1. 

The fermion determinant is real both for (3 = 1 and (3 = 4. This is obvious for (3 = 1. 
For (3 = 4 the reality follows from the identity q* = 02^2 for a quaternion real element q, 
and the invariance of a determinant under transposition. We thus conclude that Monte- 
carlo simulations can be performed for (3 = 1 and (3 = 4. For these values of (3 the 
partition function can also be interpreted in terms of quarks and conjugate quarks and 
chiral symmetry will be restored for arbitrarily, small nonzero /1, whereas a condensate of 
a quark and a conjugate anti-quark develops. Indeed, this phenomenon has been observed 
in the strong coupling limit of lattice QCD with two colors 
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Alternatively, one may write det T> + m = det 70 (T> + m) with the fermion matrix [61 




Let us study this model for m = 0. In terms of the eigenvalues, of the baryon 
number density is simply given by 



1 



n 



E 



1 1 

+ 



A fc + /i \* k + n 



(12) 



A natural electrostatic interpretation of this result is that n q is the electric field at location 
fi due to charges located at the positions of the eigenvalues. As was shown by Ginibre 



Bz] , the eigenvalues of W are distributed homogeneously in a disk in the complex plane. 
For small /1 we thus have that n q ~ fi in the quenched approximation. 

For unquenched QCD we expect that n q = for /j < fi c 7^ 0. By the electrostatic 
analog this is only possible if there are no eigenvalues in the domain /1 < \x c . The effect of 
the fermion determinant is to average the spectral density of W to zero in this region. In 
the random matrix model the phase of the fermion determinant does not result in a zero 
baryon density below /i c . Instead, n q becomes negative. The reason is that chRMT is a 
schematic model for the chiral phase transition in which the low temperature behavior 
has not been implemented correctly. The claim is that taking into account all Matsubara 



frequencies will restore this deficiency pi 



This RMT partition function at \i 7^ has inspired a great deal of interest. For recent 
literature on this topic we refer to p3| , |6^ , |5^ ,p[] . 

5.1. Yang-Lee Zeros of the Partition Function 

In this section we discuss the distribution of Yang- Lee zeros in the complex \i plane. 



In particular, we compare the zeros of the unquenched partition function for Nf = 2 and 
of the partition function with one flavor and one conjugate flavor (which we will denote by 
Nf = 1 + 1*). A recent controversial question raised in |)fj is whether the phase quenched 
partition function with Nf = 1 + 1* can teach us something about the full theory. Zeros of 
polynomials have been studied elaborately in the mathematics literature. At this moment 
we only mention a limit theorem for the distribution of zeros of polynomials with random 
coefficients which states that it converges to the complex unit circle (see for example pHU). 
Since the chRMT partition function is a polynomial in [i it can be factorized as 

Z = H(n-fi k ), (13) 

k 

where the are the complex zeros of the partition function. The baryon number density 
is thus given by 



1 ^ 1 



n. 



Q 



V E > (14) 

N I li -lik 



and can be interpreted as the electric field at the position /x from charges located at ^ 
in two dimensional complex /i-plane. 
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In QCD at T = 0, we expect that n q = for /i < fi c which is only possible if, in 
the thermodynamic limit, the number of zeros with absolute value less than /i c become 
negligible with respect to the total number of zeros. At \x c we expect a first order phase 
transition with a discontinuity in the baryon number density. In terms of the zeros, this 
implies that they should coalesce into to a cut for iV — > oo. In order to have n q = for 
|ju| < fic, this cut has necessarily to be a circle of radius fi c . Since the low-temperature 
limit has not been implemented correctly in our schematic model we will find a different 
curve for the location of the zeros. 

In terms of the cr-representation the random matrix partition function for Nf flavors 
can be rewritten 

Z N f(m = 0,n) = J Pcrexp[-n£ 2 TW t ]det n (t7t7 t - /i 2 ) (15) 

where a is an arbitrary complex Nf x Nf matrix. Since the value of the chiral condensate 
is given by the expectation value of a, the logarithm of the integrand can be interpreted 
as a Landau-Ginzburg functional for the order parameter. From a saddle-point analysis 
it follows that this partition function describes a first order phase transition ||. Since 
the equality of the real part of the free energies on opposite sides of the phase boundary 
imposes one condition on the value of fi in the complex plane, the phase boundary is 
given by a curve in the complex plane fTDfl. In the thermodynamic limit the zeros that 



are located on this curve should join into a cut. The result for this curve given by |||| 

Re^ + log^ 2 )] = -1 (16) 

is represented by the solid curve within the complex unit circle in Fig. (0). It does not 
depend on the number of flavors for n — > oo. 

The cr-model representation of the partition function for one flavor and one conjugate 
flavor is given by 

Z 1+1 * (m = 0, fi) = Jvaexpl-n^Traa^deCiak^k- /j, 2 ), (17) 

where a is an arbitrary complex 2x2 matrix and k = 03. This partition function 
is an exact rewriting of the partition function (|2|) for Nf = 2 with (pTOj) but with the 
determinant replaced by its absolute value. For even n this partition function is invariant 
under o — > au\ and // — > As a consequence it can be factorized as a polynomial in /z 4 . 

The numerical results for the zeros of the partition function for Nf = 2 (left) and 
Nf — 1 + 1* (right) are depicted by the full circles in Fig. (H). In both cases the total 
number of zeros is equal to 320. They were obtained with the help of a multi-precision 
package [[H]]. Typically, we performed our computations with about 500 significant digits. 
In both cases the complex zeros duplicate the result for the thermodynamic limit. For 
two flavors, the zeros join into two curves which are separated by ~ f° r l ar g e n 

and approach the analytical result (|16"D . Remarkably, half of the zeros coincide almost 
exactly with (|T5]). The picture is quite different for Nf = 1 + 1*. In this case we observe 
doubly degenerate zeros that approach (|16D , that approach ([IB]) rotated by 90 degrees, 
and isolated zeros that are located along the diagonals. The fraction of zeros in each 
of these three sets remains finite in the thermodynamic limit. In this figure, the total 
number in each set is given by, 120, 120 and 80, respectively. Of course, this pattern 
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Figure 2. Yang-Lee zeros of the RMT partition function with Nf = 1 + 1* (left) and 
N f = 2 (right). 

reflects that for Nf = 1 + 1* the polynomial is a polynomials in /i 4 . The vertical axis 
represents the real fi axis. For Nf = 2 we observe a second order phase transition at 
He ~ 0.52. This result can also be obtained by an extrapolation of the distribution of 
the zeros of the Nf = 1 + 1* partition function. We conclude that this partition function 
contains information of the full theory. For Nf = 1 + 1* the density of the zeros decreases 
linearly with the distance to the imaginary axis which is typical for a second order phase 
transition. 

5.2. Phase Diagram 

In the random matrix model one can include a schematic temperature dependence by 
the substitution jj, — > fi + iT. In this case the a- model representation of the partition 
function is given by |72 



Z(m,fi,T) = J doe 
where C(o) is given by 



-NC(a) 



Tr 



ln{[(cr + m)((T t +m)-(ji + iT) 2 ] ■ [(a + m)(a t + m) - (jj, - iT) 2 ]}] 



and a is an arbitrary complex Nf x Nf matrix. This model was studied for jj, — in 
P,|7|,|73| . It was found that it describes a second order phase transition at a nonzero 
temperature. Multicritical behavior was studied in [74-77]. At zero temperature, a first 
order phase transition || was found at fi ^ 0. 

Since, the value order parameter is given by the expectation value a, this partition 
function can be interpreted as a Landau- Ginzburg functional. A similar functional can be 
derived from a Nambu model It possesses a tricritical point at m = 0, \i = /i 3 and 
T = T 3 . The critical exponents at this point are given by the universal mean field critical 
exponents. We emphasize that the upper critical dimension at the tricritical point is equal 
to 3 so that up to logarithmic corrections, the mean field critical exponents are exact. 
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For a detailed discussion of the phase structure of this partition function and the physics 
of the tricritical point |79] we refer to the talk by M. Stephanov []5D| in this volume. For 



other studies of the RMT partition function at nonzero fx and T we refer to [54 



5.3. Random Matrix Triality at Nonzero Chemical Potential 

In this section, we provide evidence that the RMT Dirac operator at fi ^ shows 
universal behavior as well. In particular, we will consider the chRMT partition function 
for all three values of (3 and show that the anti-unitary symmetries of the Dirac operator 
result in characteristic eigenvalue correlations. 

Numerical simulations have been performed for all three classes. A cut along the 
imaginary axis below a cloud of eigenvalues was found in instanton liquid simulations 
Tl[] for iV c = 2 at fi ^ which corresponds to (3 — 1. In lattice QCD simulations with 

2 |S2j a depletion of eigenvalues along the imaginary axis was 

3 the eigenvalue distribution did not show any pronounced 



staggered fermions for iV c 
observed, whereas for N c 
features f|55[ . 
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1, and the imaginary parts (y) of the eigenvalues of 



In the quenched approximation, the spectral properties of the random matrix Dirac 

izing a set of matrices with 
for the eigenvalues of a few 



operator (10) can easily be studied numerically by diagonalizing a set of matrices with 
probability distribution (g). In Fig. (|3|) we show results j§3 
100 x 100 matrices for /i = 0.15 (dots). The solid curve represents the analytical result for 
the boundary of the domain of eigenvalues derived in || for (3 = 2. However, the method 
that was used can be extended [B3| to (3 = 1 and (3 = 4 and with the proper scale factors 



we find exactly the same solution. 

For (3 = 1 and (3 = 4 we observe exactly the same structure as in the previously men- 
tioned (quenched) simulations. We find an accumulation of eigenvalues on the imaginary 
axis for (3 = 1 and a depletion of eigenvalues along this axis for (3 = 4. This depletion can 
be understood as follows. For /i = all eigenvalues are doubly degenerate. This degener- 
acy is broken at fi ^ which produces the observed repulsion between the eigenvalues. 

The number of purely imaginary eigenvalues for (3 = 1 scales as y/~N and is thus not 
visible in a leading order saddle point analysis. Such a \J~N scaling is typical for the regime 
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of weak non-hermiticity first identified by Fyodorov et al. [ 84 ] . Using the supersymmetric 
method of random matrix theory the \/N dependence was obtained analytically by Efetov 



35 1 . Obviously, more work has to be done in order to arrive at a complete characterization 



of universal features [BE] in the spectrum of nonhermitean matrices. 



6. CONCLUSIONS 

There is strong evidence from lattice QCD simulations and partially quenched chiral 
perturbation theory that eigenvalue correlations below a scale of F 2 /T I \^V are given by 
chRMT. We emphasize that this is an exact result. On the other hand, qualitative 
results have been obtained for a schematic chRMT model for the chiral phase transition 
at nonzero chemical potential and temperature. We have discussed the mechanism of 
quenching and the distribution of the Yang-Lee zeros and found that the phase quenched 
partition function contains information of the critical chemical potential of the full theory. 
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